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Oh' energies in the short-range potential V{x) = ap^{x) + Ap^(x) with p(x) = sechax 



X 



Within the framework of Lanczos-inspired perturbation theory wave functions and 
energies in the short-range potential V{x) = ap^{x) + Ap^( 
and a small coupling A are shown obtainable in closed form. 



PACS 03.65.Ge 02.30.Gp 02.30.Hq 03.65.Db 



^ norbe.tex, February 7, 2008; e-mail: znojil@ujf.cas.cz 



Poschl- Teller [|l| potential V^^\x) = — + l)secli^a x resembles the celebrated 
harmonic oscillator. Within the so called shape invariant family these two spa- 
tially symmetric exactly solvable potentials form a unique subset. In perturbation 
theory there seems to emerge one of the most significant differences between them. 
In contrast to an enormous interest in the various anharmonic forces ^ there ex- 
ists virtually no analysis of a perturbed Poschl- Teller model in the current literature. 
Partially, we intend to fill the gap. This letter is concerned with the most elementary 
quartic example 

^{PPT){x) = -l^lppT)i'{PPT){x). (1) 

The parity is conserved, 'ip{ppT){.—x) = {.—^Y4^(ppt){.x), p = 0, 1, and all the un- 
perturbed A = bound states are available in closed form. Unfortunately, after we 
re-scale a — 1 for simplicity we immediately notice that the normalizable states 
form a mere finite set with /t(ppy) = ~ 2A^ ~ P = f^{N^p) and wave functions 

tanh'^ X / 1 1 \ 

(-I'^S.,) = ^^Ji^ ■ .F, - AT + -. -TV. 1 + .(7V.P). (2) 

where < 2 + p < ^ . This is the reason why these mutually orthogonal 
elementary functions cannot be used as an unperturbed basis. 

A key to our present A 7^ construction will lie in the use of a non-orthogonal 
basis. We shall employ the ansatz 

il>(ppT)ix) = tanh^x ^rx^^ (3) 

^ ^0 cosh^'^+'f^^^) X 

with K(^ppT) = K{N,p)+2e{\). Its use is inspired by the general method of Lanczos |^ 
and its perturbative implementations ||^. The basis itself is taken from the particular 
terminating expansions @ which define the unperturbed coefficients c^^-* such that 
Cy^^ 7^ and Cq^ = 1 while c^^+i = = ■ ■ ■ = 0. 

As always in similar constructions one inserts eq. (|^) in our differential eq. (|I]). 
This gives the current recurrence relations which may be written in the three-term 
form 

A c„_i + pnCn + a„+ic„+i = 0, n = 0, 1, . . . (4) 



fijfi+l) ^ 4 A 
dx'^ cosh^ a x cosh^ a x 
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or as an infinite-dimensional linear algebraic problem 
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In general, unfortunately, equation (§) cannot be treated as an infinite-dimensional 
limit of its truncated matrix subsystems due to their non-variational, power-series 
origin |]^, ||. Benefits brought by the tridiagonality concern only the wave functions. 
Their coefficients are defined by the closed determinantal formulae M 
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0,1, 



(6) 
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Once we postulate, in the spirit of the current perturbation theory. 



.(A) = c(°) + Ac« + A^c(^) + ... 



n 



1,2,... , 



(7) 



we may generate the separate coefficients from eq. by elementary algebra. This 
type of dependence of the wave functions upon the energy which is not known re- 
sembles the Brillouin-Wigner perturbation method able to treat the energy as an 
external parameter. In the similar vein the use of the variational energies has been 
recommended in the semi-analytic power-series context, e.g., by Hautot and Tater 
||10|| . The subtle problem of the determination of energies in the present perturbative 
example is to be settled in what follows. 

For the sake of definitness, perturbations with c„(A) = c^°^ + Xhn{X) will be 
normalized to /iq = 0. In order to reduce the complexity of formulae we shall only 
pay attention to the ground state problem with quantum numbers p = and = 0. 
This implies that for the sufficiently small perturbations our bound state exists at 
any non-negative /i > The details of transition to p = 1 and/or to > 
\T^ are not too interesting in the present context. For our choice of p = and 
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= the accepted normalizations simplify the very first recurrence relation. It 
reads Po + \ aihi = 0, defines the coefficient hi = 0{1) "to all orders" and suggests 
a change of the notation, /3o = A 70. As long as 70 = 0{1) we may recall our explicit 
Pq = e {e + fi + 1/2) and infer that e = e{X) = 0{X). Using a particularly convenient 
"strength-reparametrization" a = (2/i + l)/4 we shall postulate that 6{X) = r(A) — a 
with some not yet known analytic function. In terms of its expansions, say, 

r(A) = a + 6A + cA^ + dX^ + fX^ + gX^ + O [X^) (8) 

our new "energy-parameters" a,b, . . . enter the matrix elements a„ = n/2 — n'^ — 2nT, 
/3o = — and /5„ = ra^ + 2nT + /Sq. For illustration we may insert our ansatz in 
Po = 2 abX + (2 ac + 6^) + O (A^) giving 70 = 2 a6 + (2 ac + 6^) A + O (A^) etc. 

Wave functions have already been specified by eq. @. Apparently, there are 
no conditions imposed upon the energies. This is a paradox which we are going to 
explain now. In the first step we insert the coefficient hi = — 70/ai in the first, second 
and third row of recurrences (|) or (P). The former relation becomes an identity but 
the next one preserves a genuine three-component character, 

2{2e + iJ + 2)h2 = l + {e + l){e + IJ + 3/2) hi. (9) 

Up to a tiny perturbation the infinite rest with n = 1, 2, . . . has just a two-term form 

{n + 2){2e + n + fx + 2) hn+2 = {e + n + l){e + n + + 3/2) hn+i + A /i„ (10) 

giving immediately the asymptotic estimate valid for all |A| <^ 1, 

, T{e + n)T{e + n + ^1 + 1/2) ./^ ^ 

This implies the convergence of ip{ppT){x), up to the central a; = 0, rigorously. 

The boundary of the circle of convergence coincides with the centre of the spatial 
symmetry. Hence, for our particular even-parity choice of p = the first derivative 
of the wave function must have a nodal zero there, 

, , , ^ (2n + k) /i„ , , 

^xV'CPPT) (a;) ~ smhx 2^ — 2^+^+1 ^ 0, x ^ 0. (12) 

n=o cosn X 
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As long as n /i„ ~ n^^/^ the infinite sum itself is divergent in the origin. The whole 
expression ( ]T^ ) exhibits an x oo indeterminacy there. The conservation of parity 
must be imposed "by brute force" . Vice versa, the strict validity of the boundary 
condition dx4'(ppT){0) = represents precisely the "seemingly lost" quantization 
condition. 

In the nearest vicinity of x = the finite but very large value of the sum in eq. ([T2|) 
is equal to a positive real number multiplied by the coefficient /i2- This observation 
follows from eq. (p!0|) which gives, step-by-step, h^^O) = 2(/i + 5/2) /(3yU + 9) ■ /i2(0), 
h^(0) = 3(/i + 7/2)/(4/i + 16) ■ /i3(0) etc. In the interval e > -(/i + 3)/2 all the 
multiplication factors remain positive. This guarantees that the parity is conserved, 
in a semi-infinite interval of A, if and only if the coefficient /12(A) vanishes in the 
leading order approximation. 

This is one of our most important observations. The required leading-order 
change of sign of /12(A) (i.e., relation h'^'^ = 0) may be interpreted as the consequence 
of the Sturm-Liouville oscillation theorems [jl2l. In their light, for an increasing 
or decreasing leading-order energy E = —k^ = —[2a — 1/2 + (9(A)]^ the nodes of 
dx'^{ppT){x) would smoothly move along the real axis of x. In this sense the condi- 
tion h'^^ = of coincidence of one of the nodes with the origin reflects the zero-order 
physical interpretation of the flrst coefficient a in our expansion (|^) and forces us 
to postulate, consequently, the disappearance of all the subsequent leading-order 
coefficients, 

h^^ = hll^ = . . . = 0. (13) 

Before writing it down as an explicit algebra let us first move up to the next order 
approximation. Due to the danger of a re-introduction of the asymmetry (or, in the 
other words, of a spike-shaped discontinuity in ip(ppT){x) at x = [|lTl) in the higher 
orders of A we just have to repeat the previous argumentation. In the second order 
we must first notice that the vanishing condition (|TB|) changed also our recurrences 
(^). The role of the "last genuine three-term" relation (without a guarantee of a 
sign-preservation) moves from n = 1 to n = 2 in eq. (|^). As long as our estimate 
(pj]) remains valid within its own error bound l + 0{l/n) the role of an overall sign- 
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determining factor A h2{X) = (9(A^) is taken over by the new norm A h^iX)) = (9(A^). 



Mutatis mutandis we get the second-order condition h 



(2) 
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(2) 



general, 



'^n{X) = A h^X) = A"/ii") + A"+i/ii"+i) + . . . , 



k 



n = l,2,... 



and, in 
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This is a perturbative generahzation of the standard termination rules. It modifies 
our above naive expectations and incorporates correctly the physical definite-parity 
requirement in our original ansatz In a more compact notation let us put cq = 
/o(A) = 1 and A/i„(A) = A"/„(A) with /„(A) = (9(A°). This modifies our system of 
equations _ 

Po Xai 
X XPi A^a2 

A2 X^(32 A^as 



\ f f \ 

Jo 



X' A3/33 
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/i 



0. 



(15) 
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All these equations may be solved by the closed determinantal formula again. 



/n+l(A) 



T{fi+l + 2Xri) 
r{fi + n + 2 + 2Xr]) 
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(16) 



Mn{X) = det 



(3o{z) Xai{z) 
X X(3i{z) X'^a^iz) 
A2 

X^-'Pn^i A"a„ 

A" X^(3n 

It pre-factorizes the wave function coefficients in a certain less usual way. We may 
re-shuffle it slightly once more, in the more traditional Pade- approximation spirit, 
with /o = 7-1 = 1 and 



V 



7i-i 



j = 1,2,... . 



(17) 



^ (-ai)(-a2)...(-aj) ' 
This representation of the wave function coefficients is related to our final, "optimal" 



5 



/5o -A 
-«! Pi —A 

-02 P2 



. 



1^ 



recurrences 

7o 

-A 7i 

V ••■•■■•Aw 

The proof of their equivalence to our previous equations is trivial One has only to 
verify that the new notation is not inconsistent since the first row of eq. ([T8| ) just 
reproduces our old leading-order "change of the notation" Po = X 70- 

We are now prepared to compute the unknown auxiliary variables b,c, . . . and, 
in effect, the perturbed energies. Firstly we get rid of the "already known" 70 in 
eq. (p!8|). This only means that we omit the first line and abbreviate A2 = a2lo and 
Bi = /3i7o — tti in the remaining equations 

/ A -A \ 

-A2 P2 -A 

-as -A 
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Now we may parallel the above study where, retrospectively, Bq = Po- In the other 
words we only have to require that Bi = A 71 vanishes in the unperturbed limit 
(cf. the first line in eq. (|19|)). In contrast to its trivial j = predecessor the new 
j = 1 situation imposes a constraint upon Bi = 2a + 1/2 + 2ab{l + 2a) + 0{X). Its 
zero-order component must be zero. This determines the physical value of the first 



moment b = b{a) = r 



(1) 



-(1) 



1 + 4a 



(20) 

4 a (1 + 2 a) ^ ^ 

As expected on variational grounds this result means a positive first order change in 
our energy E = —k^ where k = /t(A) = [4r(A) + 2yU — l]/4. 

The whole j j + 1 procedure can be, obviously, iterated. We arrive at the 
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equations 
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with the elements A^+i = aj+i7j_i and Bj = Pj'jj-i — aj-yj^2- For any j, as we 
have aheady shown, the first-hne rule Bj = \jj = 0{\) determines the value of 
the moment t^^^ . Beyond the formal choice of j = (confirming the consistency of 
a = T^^^ = /x/2 + 1/4) and of j = 1 (giving b = t^^^) we must perform symbolic 
manipulations which become a task for the computer. Without its assistance one 
easily obtains just the next and quite compact j = 2 formula 

;i+4a)(8a2 + 8a3 + 3a+l) 



c = r 



(2) 



-1/32 



(1 + 2a) a3(l + a) 

With the assistance of the MAPLE language |jl3| the evaluation of the next correc- 
tions remains entirely straightforward giving d or 

(l + 4a)p7 



-(3) 



128 (1 + 2 af (1 + af (3 + 2 a) 



with 



I?7 = 3 + 26 a + 97 + 168 + 196 + 288 + 320 + 128 
or / or rather 

^(4) 



:i + 4a)I? 
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with 



2048 (1 + 2 af (1 + (3 + 2 af (2 + a) 



Vi2 = 90 + 1335 a + 8815 + 32715 + 69135 + 54250 



-106568 - 340152 a' - 378096 - 165184 + 22272 + 43008 a" + 10240 a 

etc. These results inspire the general ansatz 

2-M{K) (1 +4a)p^ 



12 



r 



{K) 



^2^-1(1 + 2a 



2ii'-l 



l)^"'(a + 3/2)^"' (a + 2)^"^ 



. (a + KI2) 
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re-confirmed by the explicit evaluation of the next-order correction g, 

^(5) ^ (l + 4a)Pi8 

8192 a9 (1 + 2 a)^ (1 + a)^ (3 + 2 af (2 + (5 + 2 a) 

with the rather complicated 

Vis = 3780 + 80892 a + 794817 + . . . + 458752 a^^. 

This did not disprove the obtrusive hypothesis that are polynomials in a with 
integer coefficients and with a growing degree L = L{K) = {K + 1){K + 2)/2 — 3. 

We may summarize our approach to eq. (|I]) as a scheme which can be easily gen- 
eralized. Its use of a suitable non-orthogonal basis may find applications in different 
settings. In a way paralleling the Lanczos method (with a strong numerical flavor) 
the idea itself has already been tested on several "less solvable" semi-numerical ex- 
amples [D, |ll| . In contrast to them our present construction seems to exhibit a much 
closer similarity to the current perturbation constructions of anharmonic oscillators. 
Even in this comparison, our new example seems to hide several pleasant surprises 
(e.g., a possible non-zero radius Xmax of convergence) which could make it worth of 
a further study. 
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